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Abstract 

Invariant p j spectra and rapidity densities covering a large rapidity range (-0. 1 < y < 3.5) are presented for 1& and K ± mesons from 
central Au+Au collisions at yjs^N = 62.4 GeV. The mid-rapidity yields of meson particles relative to their anti-particles are found 
to be close to unity (tt~ /n + ~ 1, K~/K + ~ 0.85) while the anti-proton to proton ratio is p/p ~ 0.49. The rapidity dependence of the 
n~ /n + ratio is consistent with a small increase towards forward rapidities while the K /K + and p/p ratios show a steep decrease 
to ~ 0.3 for kaons and 0.022 for protons at y ~ 3. It is observed that the kaon production relative to its own anti-particle as well 
as to pion production in wide rapidity and energy ranges shows an apparent universal behavior consistent with the baryo-chemical 
potential, as deduced from the p/p ratio, being the driving parameter. 

Keywords: heavy ion collisions, strangeness enhancement, baryon chemical potential 
PACS: 25.75q 
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Introduction 

As the collision energy has increased with the advent of 
new relativistic heavy-ion accelerators, from AGS energies 
( V*aw ^ 4.9 GeV) to those achieved with the SPS ( y«AW ^ 
17.3 GeV) and recently with RHIC ( V^w < 200 GeV), the 
system created in heavy-ion collisions has been found to evolve 
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from one that is baryon rich to one dominated by mesons 
id S 01- This change is evident in the increase of cen- 
tral rapidity densities of emitted pions and kaons accompanied 
by the simultaneous shift of the net-baryon peak one or two 
units of rapidity down from beam rapidity. Due to this shift, 
the net-baryonpeak gradually moves from mid-rapidity (AGS 
and SPS) % towards high rapidity (RHIC) Jl, leaving a 
relatively net-baryon poor region at mid-rapidity at the highest 
RHIC energy. 

At AGS energies, the observed ratio of strange particles to 
pions in A+A collisions is larger than that measured in either 
p+p or p+A collisions and this ratio increases significantly with 
beam energy This behavior is understood within cas- 

cade models as arising from hadronic rescatterings involving 
heavy baryon resonances 110. [Till . In the SPS energy range it 
is conjectured fl2ll that the nuclear fireball undergoes a phase 
transition from bound hadronic matter to a deconfined quark- 
gluon plasma state (QGP). This conjecture is supported by the 
observation that excitation functions of the K + A/n, and 
H~ jn ratios at mid-rapidity are found to peak around yjs^N = 
7.6 GeV H 1 3fl and then decrease slightly with increasing energy. 
The energy dependence of other observables, like the onset of 
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a plateau for the inverse slope parameter of kaon spectra and a 
kink in the 4tt pion yields normalized to the number of partic- 
ipant nucleons N„/N parl , is also used to support the transition 
to a deconfined phase at yswv = 7.6 GeV ill 311 . The gross 
features of the K/n excitation function can be described in the- 
oretical models, but additional features are needed to describe 
the data quantitatively. Phenomenologically motivated thermal 
models can explain the observed behavior with the assump- 
tion of a QGP formed in the early stages of the collision 11211 . 
However, more recent thermal models are able to reproduce the 
'horn' structure of the excitation function by including higher 
hadronic resonant states without the assumption of a phase tran- 
sition [ 14]. The K I if ratio has a monotonically increasing en- 
ergy dependence and is always lower than the K + ratio be- 
cause of the larger share of strange quarks, which form baryons 
(mainly A's) 1 1311 . From the top SPS energy to the top RHIC 
energy the K + /n + ratio at mid-rapidity is fairly constant, sug- 
gesting that the nuclear medium reaches chemical equilibrium 
with respect to the production of strange quarks [0] . 

We present invariant pj spectra of pions and kaons in cen- 
tral (0 - 10%) Au+Au collisions at ^s NN = 62.4 GeV mea- 
sured by the BRAHMS spectrometers! 15]. At this energy it 
is possible by varying the spectrometer angles to achieve an 
overlap in terms of the baryo-chemical potential (as reflected in 
the pi p ratio) of the RHIC and lower energy SPS results. The 
data cover the rapidity interval -0.1 < y < 3.5, which comes 
close to the beam rapidity {ybeam = 4.2). Total yields deduced 
from the spectra are used to calculate particle ratios that can 
be compared to the lower energy SPS results and to theoretical 
models. The large rapidity coverage allows very different nu- 
clear media to be probed: at mid-rapidity the fireball is dense 
and the yields of produced anti-hadrons and hadrons are simi- 
lar (7T/7T+ ~ 1, K-/K + ~ 0.85, pip ~ 0.45), while at forward 
rapidity, in the fragmentation region, we observe a net-baryon 
rich medium with 4 times smaller meson densities and a 20 
times smaller pip ratio. By covering a large interval for im- 
portant physical parameters {e.g. rapidity, particle density, net- 
baryon density, baryo-chemical potential), these measurements 
will provide additional constraints to those models that support 
the conjecture of a phase transition to quark-gluon plasma. 

Experimental Setup and Data Analysis 

The BRAHMS experiment consists of global event character- 
ization detectors and two spectrometer arms. Collision central- 
ity is determined using a hybrid array of Si strip detectors and 
plastic scintillator tiles located at the nominal interaction point. 
The Mid-Rapidity Spectrometer (MRS) and the Forward Spec- 
trometer (FS) are small solid angle spectrometers that can be 
rotated in the horizontal plane so as together be able to achieve 
polar angular coverage from 90° to 2.3°. For the present stud- 
ies, the MRS spectrometer was positioned at 90°, 45° and 40° 
and the FS at 6°, 4° and 3° at a few magnetic field settings, due 
to the limited running time. 

Particle identification (PID) is achieved in both spectrome- 
ters using time-of-flight walls (TOFW in MRS and H2 in FS). 
In the FS a ring imaging Cherenkov detector (RICH) is used 
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Figure 1: (Color online) Spectra for K + (left) and;r + (right) from 0-10% central 
Au+Au collisions at ^s^n = 62.4 GeV in selected rapidity slices. Each spec- 
trum is multiplied with a factor of 0.2" for better visibility. For the kaon spectra, 
the values of n are 0,1,2,3,4,5,6 for the rapidities y = 0.0,0.7,0.9,2.6,2.7,3.2 and 
3.3, respectively. For the pion spectra, the values of n are 0,1,2,3,4,5,6 for the 
rapidities y = -0.1,0.1,0.8,1.0,3.1,3.2 and 3.5, respectively. The fits are mj ex- 
ponentials for kaons and pr power laws for pions. The numbers on the plot 
indicate the rapidity corresponding to each spectrum. The error bars represent 
statistical errors. 



in addition. A detailed description of the technical capabilities 
of the experimental setup can be found in references III 611 and 
fnil . The identification of charged pions and kaons with the 
time-of-flight detectors was done by using cuts in (m 2 , p) space, 
where m 2 = p 2 (tj OF /L 2 - 1/c 2 ) is determined for a given mo- 
mentum p by using the time of flight (frof) and the track path 
length(L). Momentum dependent three-cr cuts about the mean 
m 2 of a given species were applied with the additional condition 
that the level of contamination from other species at a given m 2 
and p is limited to 5% inside the cut. Clean n - K separa- 
tion is achieved up to momenta of 1.5 GeV/c in the MRS and 
3 GeV/c in the FS. For higher momenta the 3<x curves overlap 
and the contamination condition gives asymmetrical cuts in the 
(m 2 ,p) space, necessitating the use of momentum dependent 
cuts and corrections. For the TOFW, the momentum dependent 
PID corrections ranges from at 1.5 GeV/c up to 20% for pi- 
ons and 50% for kaons at 2.2 Ge V/c. In H2, the same correction 
ranges from at 3 GeV/c to 10% and 80% at 4.5 GeV/c for pi- 
ons and kaons, respectively. In the FS, the RICH was used to 
identify pions and kaons at higher momenta. The pions are sep- 
arated from kaons by comparing the reconstructed Cherenkov 
ring radii for momenta with the calculated value for pions in 
the momentum range of 2.5 to 20 GeV/c. Kaons are identified 
by the RICH above 10 GeV/c where the RICH efficiency satu- 
rates at a value of 97% Between 4.5 GeV/c and 9 GeV/c 
an indirect method labels as kaon particles those that give no 
signal in the RICH but are also not identified as protons in H2. 
The contamination of kaons with pions unresolved in RICH and 
protons unresolved in H2 depends on momentum and on rela- 
tive particle abundances. For K the contaminant contribution 
is estimated to be almost constant at 20% whereas for K + it 
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Table 1: dN/dy values for charged pions as a function of rapidity. The given errors are statistical only. The pj range specified represents the fit range. Fiducial yield 



is the integrated yield covered by the data. 
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Table 2: dN/dy values for charged kaons as a function of rapidity. The errors are statistical only. The pr range specified represents the fit range. Fiducial yield is 



the integrated yield covered by the data. 
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varies from 14% to 40%. 

Invariant differential yields, P jdp T dy » were constructed for 
each spectrometer setting and were corrected for geometrical 
acceptance, tracking and PID efficiency, contamination, in- 
flight weak decays and multiple scattering effects by using a 
Monte Carlo calculation simulating the geometry and tracking 
of the BRAHMS detector system. The feed-down correction 
for charged pions originating from the weak decays of and 
A was applied as described in reference 0] and amounts to 5% 
of the measured yield in MRS settings and 7% in FS settings. 
The feed-down correction was applied directly to the pj inte- 
grated rapidity densities. By merging all of the spectrometer 
magnet settings, invariant pj spectra were extracted in several 
rapidity intervals and fitted with different functions in order to 
extract the integrated yields (see Fig.[TJ. We estimate the point- 
to-point systematic errors to be <5% and the systematic errors 
from normalization, tracking efficiencies and other corrections 
to be ~8%. The spectrometer acceptances allow measurements 
down to pj = 0.2 GeV/c for pions and pj = 0.25 GeV/c for 
kaons, with the absence of lower momentum data contributing 
to the systematic errors for the integral yield. This uncertainty 
is largest for pions because of their lowest average pj. In order 
to estimate the extrapolation uncertainty, we used both a power 
law distribution of the form A( 1 +prl Po) B and an mj exponen- 
tial function to fit the pion spectra. Of these two, the power law 
distribution gives the best fit to the experimental results over the 
observed momentum range, and for very low pj < 0.1 GeV/c 
agrees with the measurements made by the PHOBOS collabo- 
ration at y = 0.8 1 1811 . For the mj exponential function the fit 
range is limited to pj < 1 GeV/c. The kaon spectra are equally 
well described by an mj exponential function, Ae^~ mT ' T \ and by 



a Boltzmann distribution. All the presented data are tabulated 
and available on the BRAHMS website lfl9ri . 



Results and Discussion 

The resulting pj integrated yields, dN/dy, for pion and kaons 
are shown in Tables Q] and [2] The dN/dy distributions, taking 
the statistically weighted average of the two functional forms 
used in fitting each species, are shown in the upper panel of 
Fig. [2] The horizontal error bars indicate the width of the ra- 
pidity slices and the vertical error bars, smaller than the marker 
size, are statistical errors. The square brackets below and above 
each data point indicate the dN/dy values extracted with the 
employed functionals (see Tables [T] and |2]i and serve as an esti- 
mate for the systematic errors due to extrapolation to low trans- 
verse momentum. 

The bottom panel of Fig.|2]shows the rapidity dependence of 
the average transverse momentum (pr) for pions and kaons. 
For pions, {pj) ~ 0.41 GeV/c at y = and decreases to 
~ 0.32 GeV/c at y > 3.0 while for kaons, (pr) drops from 
~ 0.65 GeV/c at y = to ~ 0.5 GeV/c at y = 3.2. The drop 
of the averaged transverse momentum at forward rapidity is ac- 
companied by a relatively large decrease in particle densities 
and in collective models can be explained by the decrease of 
the collective radial flow velocity. 

Figure [3]shows the rapidity-dependent anti-hadron to hadron 
integrated yield ratios for pions, kaons and protons. The p/p 
ratios are obtained by using the yields presented in reference 
12511 . The baryon yields are not corrected for feed-down from 
hyperons (see [25]). The n~/n + ratio is approximately equal to 
unity over the entire rapidity range. The kaon and proton ratios 
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Figure 2: (Color online) Upper panel: dN/dy as a function of rapidity for jt* 
and S' ± from 0-10% central Au+Au collisions at ysjviv = 62.4 GeV. The er- 
ror bars are statistical. The square brackets below and above each data point 
indicate the dN/dy values extracted with the selected functionals (see Tables^ 
and [2] and text). The experimental uncertainties on the extrapolated yields are 
below 5% and not shown in the figure, except for the pion yields at y = 3.5 
where these amount to ~ 30% and are indicated by the gray box. Bottom panel: 
(pr) dependence on y. The statistical error bars are covered by the symbols. 



at mid-rapidity {K~/K + ~ 0.85, p/p ~ 0.49) are lower than 
the corresponding ones measured at the top RHIC energy 12411 . 
but they are still characteristic of a high degree of anti-matter to 
matter equilibration. At forward rapidity we observe a decrease 
of the kaon ratio to a value of ~ 0.35 at y = 3.3. Possible ex- 
planations include the competition between A baryons and K 
mesons for the available strange quarks and associated produc- 
tion (e.g., p + p — > p + A + K + ) which increases the number 
of positive kaons. Both of these mechanisms depend on the 
net-baryon content and, consequently, lead to a decrease of the 
K /K + ratio at forward rapidity. The p/p ratio decreases sig- 
nificantly with rapidity, reaching p/p — 0.077 ± 0.004(sfaf.) at 
y=2.3 and p/p = 0.022 + O.OOl(stat) aty = 3. 
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Figure 4: (Color online) Rapidity dependence of the K/n ratios in 0-10% cen- 
tral Au+Au collisions at ^s^n = 62.4 GeV. The error bars are statistical errors 
and the square brackets show the systematic uncertainties due to the yield ex- 
trapolation at low pr. 
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Figure 3: (Color online ) Anti-particle to particle ratios as a function of rapidity 
in 0-10% central Au+Au collisions at yi/vA/ = 62.4 GeV. The error bars are 
statistical only. 



In Fig. [4] we show the rapidity dependence of the K/n ra- 
tio. Because the rapidity intervals where the yields of the two 
species were extracted are not the same at forward rapidity, we 
used a linear interpolation procedure between the closest cov- 
ered points to obtain the meson yields for additional points in 
rapidity. We checked this procedure by assuming Gaussian ra- 
pidity distributions and found very similar results. The K + /n + 
ratio was found to be 0.159 + 0.011 at mid-rapidity and is al- 
most constant as a function of rapidity. The K /n~ ratio has a 
value of 0.13 + 0.01 at mid-rapidity and shows a steep decrease 
for y > 2.5 with a value of ~ 0.05 at y = 3.2. The different 
rapidity dependence of the positive and negative charge K/n 
ratios is similar to that found in central Au+Au collisions at 
■\/snn = 200 GeV J3l but the difference between the two ratios 
is three times larger at y=3. 

In the following we use calculations with two micro- 
scopic transport models for comparison with our data. Ultra- 
Relativistic Quantum Molecular Dynamics (UrQMD) S 
is the extension of RQMD ifiol [Till which was developed to 
describe physics at AGS and SPS energies. At low energies, 
ys/wv < 5 GeV, it describes the nuclear collisions in terms 
of interactions between known hadrons and their resonances 
while at higher energies the dominant mechanism is through the 
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Figure 5: (Color online) K~ /K + ratio dependence on the p/p ratio. The solid 
circles are from 0-10% central Au+ Au collisions at -Jsr^tj =62.4 GeV obtained 
in the present work. The open symbols are BRAHMS data from ref. l24ll and 
lower energy data from 131 123. 130t 131L l32l [3311 . The error bars represent sta- 
tistical and systematic errors. The curves are calculations with UrQMD (black 
lines) and AMPT (gray lines) for central Au+Au collisions at V.vaia>=200 GeV 
(dashed lines) and 62.4 GeV (solid lines) and for central Pb+Pb collisions at 
-y/sjv/v = 17.3 GeV (dot-dashed lines). The curve for the UrQMD calculation at 
17.3 GeV is very close to that at 62.4 GeV making it less visible. 
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Figure 6: (Color online) K/n ratios as functions of the p/p ratio. The circles 
are data from 0-10% central Au+Au collisions at -^s^n = 62.4 GeV and the 
triangles are results obtained at mid-rapidity in the SPS experiments at different 
energies flUl . The curves are UrQMD (thick lines) and AMPT (thin lines) 
calculations for central Au+Au collisions at ysjvjv = 17.3,62.4 and 200 GeV. 
The black curves are for the K + /n + ratio and gray curves are for the K~ fiC 
ratio. 



color string excitation. AMPT (A Multiphase Transport Model) 
I 22l 23] is a microscopic model developed for nucleus-nucleus 
collisions at RHIC and LHC energies. AMPT uses the HI- 
JING model for the initial space-time configuration of the par- 
tons and strings. Both UrQMD and AMPT include a Lund-type 
string fragmentation procedure followed by hadronic rescatter- 
ings. The models were run in minimum bias mode with a 0- 
10% centrality cut placed on the distribution of particles falling 
in the acceptance of our multiplicity detector (I77I < 2.2), closely 
mimicking the experimental conditions. 

Figure shows the dependence of the K~/K + ratio on the 
p/p ratio in central nucleus-nucleus collisions at energies rang- 
ing from yfsNN — 5 to 200 GeV. The data points obtained in 
this work and from the other RHIC energies are obtained in 
different rapidity slices, while the SPS and AGS points are ob- 
tained at mid-rapidity. Comparing the low and high energy re- 
sults, we observe a systematic dependence of the K~/K + ratio 
on the p/p ratio. The calculations made with UrQMD(thick 
lines) and AMPT(thin lines) for central nucleus-nucleus colli- 
sions at ysMv = 200, 62.4 and 17.3 GeV do not reproduce 
quantitatively the dependence of the K / K + ratio on the p/p 
ratio which seems to be universal over a large energy range. 
However, UrQMD calculations at the three energies give sim- 
ilar K~/K + ratios for a given p/p value. This is most likely 
due to the thermalization of the system reached via secondary 
rescatterings, which includes formation, decay and regeneration 
of many resonances [34]. 

Figure [6] shows the dependence of the Kjn ratios on the 
pip ratio. Our data points are obtained in different rapidity 
slices at the same energy, y/s^N = 62.4 GeV, while the SPS 
points are obtained at mid-rapidity in central Pb+Pb collisions 
at ys^ = 6.3,7.6,8.8, 12.3 and 17.3 GeV. In Au+Au colli- 
sions at 62.4 GeV, the fragmentation peak is estimated to be in 
the interval 2.5 < y < 3.3 12511 with a net-proton density of 
~ 30 and p/p ratio values at forward rapidity, which are in the 
same range as those measured at mid-rapidity at s/s^n = 12.3 
and 17.3 GeV |6J]. As in the case of K /K + ratio, we observe a 
smooth dependence of the Kjn ratios on the p/p ratio. The 
curves show calculations for central Au+Au collisions from 
UrQMD (thick lines) at = 17.3,62.4 and 200 GeV and 

AMPT (thin lines) at ysMV = 62.4 GeV. None of the models re- 
produce quantitatively the dependence of the K/n ratios on the 
p/p ratio but qualitatively we observe that UrQMD reproduces 
the trend of the data. 

Figure [7] displays the inverse slope parameters for kaons as 
a function of the p/ p ratio. The inverse slope for the highest 
ratio (0.45) represents the result from a simultaneous fit to the 
spectra in < y < I. Also shown are the mid-rapidity results 
from SPS. The inverse slopes of the spectra for positive and 
negative kaons measured at forward rapidity in our dataset are 
~ 20 MeV smaller than the ones measured at mid-rapidity in 
SPS experiments at the same p/p ratio. Even though the val- 
ues are consistent within error bars, the difference suggests that 
the transverse momenta are not governed by the baryo-chemical 
potential. The mid-rapidity particle densities for negative pions 
at SPS at yiiviv = 17.3 and 12.3 GeV are approximately two 
times larger than the negative pion densities at the forward 
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Figure 7: (Color online) Inverse slopes (7" c (j) for kaons as a function of the 
p/p ratio. The BRAHMS points are from 0-10% central Au+Au collisions at 
tJsnn = 62.4 GeV (this analysis) at different rapidities while the SPS points 
are from mid-rapidity at different energies I3I l3lh . The error bars represent 
statistical errors, while the square brackets show systematic errors. 



rapidities at y/s^N =62.4 GeV. Smaller radial flow from the less 
dense local system at 62.4 GeV could be the cause of the small 
difference between inverse slopes. The inverse slopes obtained 
at forward rapidity are higher for positive kaons than for the 
negative kaons. This is also observed at mid-rapidity at SPS 
energies. 

In a chemical analysis, the pj p ratio has an approximate cor- 
respondence with the baryo-chemical potential through the for- 
mula p/p — exp(-2ps/T) for a given freeze-out temperature 
T. Hence if T is the same in the two cases, this would im- 
ply that the local system formed at high rapidity at RHIC (62.4 
GeV) is chemically equivalent with the system formed at the 
two highest SPS energies at mid-rapidity, both being controlled 
by the baryo-chemical potential. The strangeness chemical po- 
tential is fixed by the baryo-chemical potential and chemical 
freeze-out temperature provided that the local net strangeness 
vanishes. A rapidity dependent baryo-chemical potential pg has 
been suggested in thermal models Ir26l 27 . 28|. The observed 



dependence of the chemical freeze-out composition in differ- 
ent colliding systems and different rapidities only on the baryo- 
chemical potential and temperature supports the idea of local 
chemical equilibration and the existence of a universal (T - Pb) 
freeze-out line. 



Summary 

In summary, we have measured the transverse momentum 
spectra and inclusive invariant yields of charged pions and 
kaons in central Au+Au reactions at 62.4 GeV. The anti- 
particle/particle ratios for kaons and protons show a steep de- 
crease at forward rapidity. The charge dependence of the K/jt 
ratio at forward rapidities is understood in the framework of 
microscopic models (i.e. UrQMD) as resulting from the as- 
sociated production in a baryon rich medium. The production 



mechanisms enhance the fraction of s quarks ending up in hy- 
perons, thus depleting the K yield. We observe in central 
nucleus-nucleus collisions a common dependence of the par- 
ticle ratios (K/n, K /K + ) and kaon spectra inverse slopes on 
pip, which reflects the baryo-chemical potential, whether mea- 
sured for different energies at mid-rapidity at SPS, or at differ- 
ent rapidities at yjs^N = 62.4 GeV. This is consistent with a 
picture where in nucleus-nucleus collisions, at a given energy, 
the local fireballs formed at different rapidities freeze out on a 
T - pb line which coincides with the T - ps freeze-out line 
previously observed in a wide energy range but only at mid- 
rapidity. The baryo-chemical potential interval covered at for- 
ward rapidity extends to high values and is equivalent to the 
two highest SPS energies, but not quite high enough to probe 
the horn structure in the K + /n + ratio excitation function. In the 
calculations made with UrQMD and AMPT models, K/n ra- 
tios are not well reproduced for large baryo-chemical potential 
found either at high rapidity or mid-rapidity at lower energies 
although UrQMD seems to reproduce qualitatively the trend of 
the data. Also, the universal dependence of the K~/K + on the 
p/p ratio is not well explained quantitatively by these micro- 
scopic transport models. However the UrQMD model seems to 
reproduce qualitatively very well this dependence. 

Acknowledgements 

This work was supported by the Division of Nuclear Physics 
of the Office of Science of the U.S. Department of Energy un- 
der contracts DE-AC02-98-CH 10886, DE-FG03-93-ER40773, 
DE-FG03-96-ER40981, and DE-FG02-99-ER41 121, the Dan- 
ish Natural Science Research Council, the Research Council 
of Norway, the Polish Ministry of Science and Higher Edu- 
cation (Contract no 1248/B/H03/2009/36), and the Romanian 
Ministry of Education and Research (5003/1999, 6077/2000). 
We thank the staff of the Collider- Accelerator Division at BNL 
and the RHIC Computing Facility for their support to the ex- 
periment. 

[1] L. Ahle et al, E802 Coll., Phys.Rev. C57 (1998) 466 
[2] L. Ahle et al, E802 Coll., Phys.Rev. CS8 (1998) 3523 
[3] S.V. Afanasiev et al, NA49 Coll., Phys.Rev. C66 (2002) 054902 
[4] I.G. Bearden et al, BRAHMS Coll., Phys.Rev.Lett. 94 (2005) 162301 
[5] L. Ahle et al, E802 Coll., Phys.Rev. CS9 (1999) 2173 
[6] C. Alt et al, NA49 Coll., Phys.Rev. C73 (2006) 044910 
[7] I.G. Bearden et al, BRAHMS Coll., Phys.Rev.Lett. 93 (2004) 102301 
[8] L. Ahle et al, E866 and E917 Coll., NuclPhys. A638 (1998) 57C 
[9] S. Ahmad et al, E891 Coll., Phys.Lett. B382 (1996) 35 
[10] H. Sorge, H. Stocker and W. Greiner, Ann.Phys. (NY) 192 (1989) 266 
[11] H. Sorge, H. Stocker and W. Greiner, NuclPhys. A498 (1989) 567 
[12] M. Gazdzicki and M. Gorenstein, Acta Phys.Polon. B30 (1999) 2705 
[13] M. Mitrovski et al, NA49 Coll., J.Phys. G32 (2006) S43 
[14] A. Andronic, P. Braun-Munzinger, J. Stachel, Phys.Lett. B678 (2009) 516 
[15] I. Arsene et al, BRAHMS Coll., NuclPhys. A7S7 (2005) 1 
[16] M. Adamczyk et al, BRAHMS Coll., Nucl.Instr.Meth. A499 (2003) 437 
[17] R. Debbe et al, NuclInstr.Meth. A570 (2007) 216 
[18] B.B. Back et al, PHOBOS Coll., Phys.Rev. C77 (2007) 024910 
[19] http://www.rhic.bnl.gov/brahms/WWW/publications.html 
[20] S.A. Bass et al, Prog.Part.Nucl.Phys. 41 (1998) 225 
[21] M. Bleicher et al, J.Phys. G25 (1999) 1859 
[22] B. Zhang et al, Phys.Rev C61 (2000) 067901 
[23] Z.-w. Lin et al, NuclPhys. A698 (2002) 375 
[24] I.G. Bearden et al, BRAHMS Coll., Phys.Rev.Lett. 90 (2003) 10 



6 



[25] I.C. Arsene et aL, BRAHMS Coll., Phys.Lett. B677 (2009) 267 

[26] F. Becattini, J.Cleymans, J.Phys. G34 (2007) S959 

[27] B. Biedron, W. Broniowski, J.Phys. G35 (2008) 044011; Phys.Rev. C76 

(2007) 054905 

[28] L. A. Stiles, M. Murray, arZ/v Jnucl-ex/ 060 1 039 A 

[29] L. Ahle et aL, E866 Coll., Phys.Rev.Lett. 81 (1998) 2650 

[30] M.van Leeuwen et aL, NA49 Coll., Nucl.Phys. A715 (2003) 161 

[31] C. Alt et aL, NA49 Coll., Phys.Rev. C77 (2008) 024903 

[32] C. Alt et aL, NA49 Coll., Phys.Rev. C73 (2006) 044910 

[33] I.G. Bearden et aL, NA44 Coll., Phys.Rev. C66 (2002) 044907 

[34] L. V. Bravina et aL, Phys.Rev. C78 (2008) 014907 



7 



